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Abstract: This paper derives explicit results for the infinite Gibbs par- 
titions generated by the jumps of an a— stable subordinator, derived in 
Pitman ||39| : |40|) . We first show that for general a the conditional EPPF 
can be represented as ratios of Fox i?-functions, and in the case of rational 
a, Meijer G-functions. This extends results for the known case of a = 1/2, 
which can be expressed in terms of Hermite functions, hence answering an 
open question. Furthermore the results show that the resulting uncondi- 
tional EPPF's can be expressed in terms of H- and G-transforms indexed 
by a function h. Hence when h is itself a H- or G- function the EPPF is also 
an H- or G-function. An implication, in the case of rational a, is that one 
can compute explicitly thousands of EPPF's derived from possibly exotic 
special functions. This would also apply to all a except that computations 
for general Fox H-functions are not yet available. However, moving away 
from special functions, we demonstrate how results from probability theory 
may be used to obtain calculations. We show that a forward recursion can 
be applied that only requires calculation of the simplest components. Ad- 
ditionally we identify general classes of EPPF's where explicit calculations 
can be carried out using distribution theory. Specifically what we call the 
Lamperti class and Beta-Gamma class. As a special application, we use the 
latter class to obtain EPPF's based on mixing distributions derived from 
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the laws of of ranked functionals of self similar Markovian excursions. The 
work serves importantly a dual purpose. One is to obtain explicit calcu- 
lations for large classes of EPPF's and hence of use to the growing num- 
ber of applications involving combinatorial stochastic processes. The other, 
perhaps surprising, is the introduction of new techniques for calculating 
explicitly certain Fox //-transforms and related quantities via probabilistic 
arguments. 

AMS 2000 subject classifications: Primary 62G05; secondary 62F15. 
Keywords and plirases: beta gamma algebra, Brownian and Bessel pro- 
cesses. Fox H- and Meijer G-functions, Lamperti type laws, stable Poisson 
Kingman Gibbs partitions.. 



1. Introduction 

Let Sa, for < a < 1, denote a positive a-stable random variable, whose law is 
specified by the Laplace transform, 

E[e"^^°] =e-^° 

for A > 0, and with density denoted as /„. Then, following (3^; 42) it is well- 
known that 

C30 

Sa — ^ ^ Ji 

i=l 

where Ji > J2 > . . . > are the ranked jump sizes of a stable subordinator 

00 

Sa{s) = JiI{U^ < s), < s < 1, with S'„(l) = Sa, 
1=1 

where the (Ui) are independent random times distributed uniformly on (0, 1) 
and also independent of the ( Ji). Furthermore, the subordinator is characterized 
by its Levy density, 

Pais) 



r(l - a) 



Following Kingman (|2lh and Perman-Pitman-Yor (|4ll : l35l : |42I ). Pitman (13£ 
discussed the laws of the ranked jumps normalized by their random total mass 
Sa and further conditioned on Sa —t, that is, 

£(Pl,P2,... I Sa^t) 

where (Pi) — [Ji/ Sa)- The distribution is referred to as the (conditional) Pois- 
son Kingman distribution with Levy density pa and denoted as VK{pa\t). Fur- 
thermore, one can create an infinite number of laws from this construction by 
mixing over t with any distribution 7 on (0, 00). The law of the sequence {Pi) is 
then referred to as the Poisson Kingman distribution with Levy density pa and 
mixing distribution 7, denoted as 

/•oo 

PK(Pa,7) = / PK(pa|t)7(dt). 
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In this paper we provide explicit calculations and interpretations for the ex- 
changeable partition probability function (EPPF) which characterizes the law of 
the exchangeable random partitions Hoo = (Hn) on N generated by the ranked 
jumps of an a stable subordinator, that is, the pk( pa,"/) partitions. Equiva- 
lently, using Kingman's paintbox representation (jSSi : llj llOt ITU ) this class of 
exchangeable random partitions can all be constructed by a random closed set 
Z C [0, 1] where Z is the scaled range of an a— stable subordinator conditioned 
on its value at a fixed time. Conditional on t this construction produces the 
PK{pa\t) partition, where the conditional EPPF was derived by Pitman (fsol ). 

Our work may be divided into two parts. One is the use of special functions 
and the theory of fractional calculus to help interpret these EPPF's and in many 
cases to obtain explicit numerical calculations. The other is to use some inter- 
esting probability distribution theory to obtain explicit results for large classes 
of EPPF's and which in turn yields calculations for various special functions. 

Specifically we will show that for general a these EPPF's may be repre- 
sented in terms of Fox iJ-functions, and for the case where a takes on ratio- 
nal values in terms of Meijer G-functions. There are several significant impli- 
cations of these representations. One is that a large number of special func- 
tions commonly appearing in, for instance, physics, probability, finance or frac- 
tional calculus, can be represented in terms of H- and G-functions. See, for in- 
stance, (jisl : [20I: [2^: 27 : 261 : 28 : 45). In addition, because these functions are well 



understood this offers additional interpretability of the relevant EPPF's. That is 
one is not merely applying a numerical calculation. The case of rational a is par- 
ticularly interesting as calculations involving general Meijer G-functions are at 
the heart of mathematical computer packages such as Mathematica and Maple. 
This literally allows one to explicitly calculate thousands of EPPF's, while in 
the present literature only a few cases of explicit EPPF's are known. However, 
as of yet, while in the general case of a we can express many EPPF's in terms of 
//-transforms there does not exist general mathematical packages to compute 
them. This brings up our other approach which is based on probability distri- 
bution theory associated with beta, gamma and stable random variables relying 
in large part on a recent work of James (jlSl ). see also (JITI)), and which further 
relies on some results in Perman-Pitman-Yor (41; 35; 42h and Pitman (|4oh. and 
some lesser known results for Sa- In this regard, the fact that the relevant com- 
ponents in the EPPF satisfy a forward recursion, which follows from a backward 
equation as seen in Gnedin and Pitman (fl^ . plays an important role. As will 
be discussed in Section [H the recursion shows that one need only calculate the 
simplest components to calculate all components via a recursion. Specifically, it 
suffices to compute the probability of having one block in a partition of integers 
{1, . . . , n}, that is, 

P({l,...,n}) = K,i[l-a]„_i, 

where is the quantity that needs to be computed, and the other notation 
will be explained shortly. We should note that even Vns was thought not to 
be easily calculated, however we will demonstrate that this can be done quite 
readily. We will also, in many cases, be able to give explicit expressions for 
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the more complex quantities. This also sets up some interesting relationships 
between various special functions, integral transforms, and probabilities. 

We provide a discussion and definition of Fox H- and Meijer G-functions 
in the appendix, which is obtained from various sources. We now proceed to 
address some preliminaries and present a more specific outline. 



1.1. Preliminaries 



Again from Pitman (|39r). for n„ an exchangeable partition of {1, 2, . . . , n} and 
a particular partition {Ai, . . . , A^} of {1, 2, . . . , n} with \ Ai\ = rii ior 1 < i < k, 
where rij > 1 and Yli=i = then letting, 



[a]n = a{a + l)(a + 2) • • • (a + n - 1) 



r(a + n) 



the conditional EPPF, associated with the PK{pa\t) partition, is defined as 
Pa{ni, . . .,nk\t) = P(n„ = {Ai, . . .,Ak}\t), where 



a\nj-l 



and one can write 



i"r(n - ka)fa{t) 
T{n - ka)fa{t) 



fa{t - v)v 



dv 



fa{tu){l — uY'' ^du 



(1.1) 



(1.2) 



Using the terminology in (|12l). call the PK{pa\t) partitions the (Q;|i)-partitions. 
Now suppressing dependence on a and 7, for each n and k, set 



Gt'Hthidt). 



Pitman (j39l ) shows that the EPPF of the pk(/9q.,7) partition is given by 



(1.3) 



Note that by setting 7 to be point mass at (|1.3p equates with (|l.ip . The most 
well-known member of this class is the case where for 9 > —a, 7 corresponds to 
the distribution of the random variable Sa,e having density 



r(6l/a + 1 



(1.4) 
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and satisfies ioi S + 9 > —a, 

Ff9-^i = + ^ r(M + i) r(g + i) 

^ "'"^ r(0/a + i) ^ " ^ r(0 + <5 + i)r(0/a + i)" 

Note Stj^o = Sa- Furthermore, the random variables satisfy the remarkable 
identity that may be fomid in Pitman (|40l . section 4.2) and Perman, Pitman 
and Yor (35). That is, for any 6 > —a, 

1 d Pe+a.l-a /, 

-FT- > (1-5) 



Sa 



->a,e+a 



where Pg+a,i~a is a beta(6'+a, 1 — a) random variable independent of Sn.fi. Note 
that the identity p.5|) plays an important role in the recent work of and 
hence, through that work, will play a prominent role here. That is, using ()1.4|) as 
the mixing density gives the EPPF of the two parameter (a, 9) Poisson-Dirichlet 
distribution, say, PB{a,9), given by 

PaAm, = ^hf,,'^^'"^ nil - (1-6) 



The quantity in (|1.6p extends beyond the stable case, as it is also defined for 
the cases of a = and 6 > corresponding to the famous result related to 
the Dirichlet process, otherwise known as the Chinese restaurant process, or 
Ewens (0, 6')-partitions. The other possibility is the case of — oo < a < and 
9 = m\a\, for m = 1,2..., referred to as (a, |a|rn) partitions. The PD(a,0) 
plays an important role in a variety of diverse applications. See Pitman ([401 ) for 
a general overview and set of references, _and in particular, its relation to Bessel 
and Brownian phenomena. See Bertoin (2) for its role in terms of fragmentation 
and coagulation phenomena. See Ishwaran and James dil [13) and Pitman (jH) 
for some applications in Bayesian statistics. 



A remarkable fact, see Gnedin and Pitman (12, Theorem 12) and Pitman (|4C 
Theorem 4.6, p. 86), is that the EPPF's generated by (frT|) . that is (fOI . and 
mixtures of the Ewens (0, 0)-partitions and (a, |a|m) partitions, constitute the 
only infinite EPPF's having Gibbs form, that is, infinite EPPF's of the form 

k 

Cn,k Y\_ ^"3 ■ 
3 = i 

This, as discussed in regards to the pd(q;, 9) family, has potential implica- 
tions both from a practical and theoretical point of view in a variety of dis- 
ciplines. In particular, the (a|t)-partitions constitute the largest and most di- 
verse of such classes. However, there are only a few examples where Vn,k has 
been computed. Besides the PT>{a,9) case, there are also the models formed by 
taking 7 as a density proportional to e~''*fa{t). In addition. Pitman ('39', sec- 
tion 8) and Pitman (|4Ci section 4.5, p. 90) show that conditional EPPF of the 
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(l/2|t ^)-partition, corresponding to the Brownian excursion partition, is such 
that G[y^\t) can be expressed in terms of Hermite functions (25, section 10.2). 
This expUcit resuh is due in part to the fact that S1/2 is equivalent in distribu- 
tion to an inverse gamma distribution with shape 1/2, and hence in contrast to 
the case of general Sa has a simple explicit density. However, clearly, given the 
fact that 7 may be quite arbitrary, it is not enough to simply know the explicit 
form of the density of /„. 



1.2. Goals and Outline 

Faced with this our goal becomes quite clear. Find methods to explicitly calcu- 
late and hopefully provide further interpretation of the quantities and 
Vn,k for general a. We wish to emphasize that we are not interested in merely 
suggesting crude numerical methods which do not have interpretability. 

Our first task will be to provide an answer to a question posed by Pitman ((40l 
Problem 4.3.3, p. 87), which goes beyond merely wanting a numerical calcula- 
tion. We paraphrase it as follows. 

Pitman l|4Cll . section 4.5, eq. (4.59) and (4.67)) shows that in the case of a = 1/2, 
the integral Ga"''°'(t) can be simply expressed in terms of an entire function of a 
complex variable, the Hermite function, which has been extensively studied. It is 
natural to ask whether Ga"''^'(i) might be similarly represented in terms of some 
entire functions with a parameter a, which reduces to the Hermite function for 
a = 1/2. 

In sections [2] and [3] of this paper we provide an answer to this question by 
showing that G^"''''' (t) can be expressed as the ratio of Fox i?-functions in the 
case of general a and Meijer G-functions in the case where a is rational. In par- 
ticular, in section [3TT1 we recover the case of the Hermite function based on the 
calculus of Meijer G-functions when a = 1/2, and show that for general rational 
a these may be expressed as ratios of sums of generalized hypergeometric func- 
tions. Additionally, in those 2 sections we show that Gi"''°''(t) can be expressed 
in terms of densities derived from Sa,ka and corresponding beta random vari- 
ables. In section |4] we obtain results for the unconditional Gibbs models, that 
is, calculations for Vn,k- In sectioi flTT] we show that one may use the calculus 
of Meijer G-functions to express many Vn,k in terms of G-functions which are 
then readily computable. Sections 14.1.11 and 14.1.21 provide new specific exam- 
ples of EPPF's. Section [5] represents our first real departure from calculations 
of EPPF's based on the theory of Meijer G-functions. In particular, we demon- 
strate that one can obtain calculations for all values of a using random variables 



connected to Lamperti (|24D . which have been recently studied in James (|l8l). 
Interesting special examples, which can be computed by various other means, 
are given in sections [5.11 and [5T2l Section[6l albeit short, is a pivotal section as it 
demonstrates the important role of recursion formulae. This, as mentioned ear- 
lier, shows that one only need to focus on calculation of the Vn,i terms in order 
to obtain the general Vn,k- As a side note, we believe that even in the case where 
Vn,k niay be represented in terms of G-functions and therefore computable, it 
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might be more efficient to use the recursion when one is interested in problems 
potentiaUy involving the calculation of all Vn,k for fc = 1, . . . , n and n — 1,2, ... . 
Such problems occur, for instance, in the implementation of generalized Chinese 
restaurant schemes as can be seen in Ishwaran and James (flih (see also Pit- 
man (fiol . section 3.1). In those cases n represents sample size of data and can 
be in the thousands. In addition the recursion can be used to obtain new recur- 
sive relationships for various quantities, we will demonstrate that in section [51 
In section [71 we describe purely probabilistic methods to calculate 14, i and in 
fact obtain expression for general Vn.k- We point out again that there are very 
few explicit examples of EPPF's so it is rather striking that we will now show 
how to obtain many of them. Section 17.11 introduces what we call the Lam- 
perti class which can be seen as an extension of section [5l Section [7.21 introduces 
what we call the Beta-Gamma class. Remarkably this class shows that Vn,k may 
be expressed in terms of expectations involving only beta and gamma random 
variables. Section [8l represents a non-trivial application of Proposition 7.4 which 
exploits the quite broad results in Pitman and Yor (Si). 

Remark 1.1. Throughout we will write Gs to represent a gamma{5, 1) random 
variable and Pa,b to represent a beta{a, b) random variable. Furthermore, un- 
less otherwise stated, we will assume that when we write products of random 
variables, that means the individual random variables are independent. 

2. Conditional EPPF 



From Schneider (fiil ) (see also (27; H)), one may represent the density of Sa in 
terms of an i7- function as follows. 



Ut) 



1 



-H. 



2,2 



(0,1), (0,1) 

Applying ()9.7|) in the appendix, one can write 

(1--,-) 



t > 0. 



(2.1) 



1 



a 



-m 



(0,1) 



t > 0, 



(2.2) 



and use this to describe Gi"'*^'' (i). 
Theorem 2.1. The function 

Gi"''=)(t) 

defined by U.^) appearing in il.l]) can be expressed as follows. 

(i) G^a'^^t) is representable in terms of ratios of probability densities as. 



r{n,k)t,-. _ 'r(fc) /a,(„,fc)(t) 

r(") f»{t) 
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where fa,{n,k) denotes the density of the random variables 

Pka,n-ka P(^k-l)a+l,7i-l-{k-l)a 

for fc = 1, . . . , n. In particular, for k ~ 1, 



(2.3) 



(2.4) 



(ii) For all Q < a < 1, G^a'^\t) is expressible as the ratio of Fox H -functions, 

(1- 1 _fc,i) - 



ttO,1 



(-n,l) 



H 



1,1 



(1--,-) 
(0,1) 



Proof. Let us proceed by first deriving the density of Sa,ka/ l3ka,n-ka- First 
for clarity we use the fact that Pka,n-ka has density, 



r(n) 



r{ka)T{n — ka) 



-ka — 1 



for < li < 1. Now setting — ka in (|1.4p . we see that the usual operations to 
obtain the density involves 



— ka-\-lj. — k 



<-'="]/„(iu)u'="-i(l-u) 
Hence due to the above cancelations, the density is 



\n—ka—l 



fa,(n,k){t) — 



Now it follows that, 



T{k)T(ji — ka) Jq 



i""" / /a(M(l-«) 



l^ka^n—ka — Pka-{'l~a.n~ka — l-}-aPka.l — a. 

Furthermore, setting = ka in (|1.5|) . one gets 



ka,l — c 



These two points yield the equivalence in (|2.3p and (|2.4p . 
In order to establish statement (ii), first write 

r(n - ka) Jo 



-ka — l 



du. 
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Substituting the expression for fa{ut) with (|2.ip . and then using (|9.10p in the 
appendix to obtain an expression for {tu)^^°' fa{ut), one sees that 



aT{n — ka) Jq 



n—ka—l 



xH. 



1,1 



tu 



(l-i-fc,i),(-fc«,l) 

{—ka, 1), {—ka, 1) 



du 



1 TjOA 



(1- i-Zc, i) 
{-n,\) 



(2.5) 



where the last equality follows from identity 2 in (j45l . p. 355) and some manip- 
ulations based on (|9.1ip . (|9.6p and (|9.7p . The result is concluded by applying 
the representation (12. 2p . 

Statement (i) of Theorem I 2 . 1 I provides a probabilistic interpretation of Gi"''^''(t). 
We shall see that, in particular, the distributional identity (|2.4p will play a cru- 
cial role in applying probabilistic arguments to calculations based on Gi"''°''(i). 
Statement (ii) expresses Gi"''''*(t) in terms of Fox i/-functions, which, among 
other things, allows one to make precise interpretations of it. Overall, Theo- 
rem l2.1l sets up a myriad of duality relationships between probabilistic quantities 
based on stable and beta random variables and a very large class of special func- 
tions. In particular, as we will show, one can use statement (i) to obtain explicit 
calculations for various Fox _ff-functions which are not yet readily computable 
by other means. Now when a = m/r, the next result shows that expressions 
in statement (ii) reduces to ratios of Meijer G-functions. The significance being 
there is that calculations of Meijer G-functions are readily available in Mathe- 
matica and other mathematical softwares. 



3. Rational Values, Products of Beta and Gamma Random variables 
and Meijer G-Functions 

When a = m/r for integers 1 < m < r, the stable random variable Sm/r can 
be represented in terms of independent beta and gamma random variables as 

follows. ( — — ) = r^Zjn,r, whcrc 



This result may be found in Chaumont and Yor (4, p. 113) (see also (jisl , section 
6)). Now, by Theorem 9 in Springer and Thompson (48, p. 733), the density of 
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Zm.,r can be represented in terms of Meijer G-functions as follows, 

( i 1 \m — 1 ~ 



where 



r-1 



^"V- = n TTir 11 



(^-i)r 

1 



(3.1) 



(3.2) 



i— 1 V ^' / j—m 

Hence the density of S22L in s is given by 

m+l 



/2 



0,r-l 



r 



.(m-|-l)/m 



(1 
(1- 



■i \r— 1 



1 j_ \m— 1 

m m 1 



(3.3) 



by absorbing the term (jn"^r~^ s~"^y"'~^'^" using (|9.10p and then applying ([97 
Note one could have used the result of ((i^ to obtain (|3.3p . However that result 
does not equate S'^/r with beta and gamma random variables. 
Now define the vectors, 



1 i + k 



i-l k 



and 



1 j \ ™ ^ / j - 1 - n 



m m r 



(3.4) 



(3.5) 



Theorem 3.1. Let m,r denote integers such that I < m < r. Then, for a 
m/r, hi. 2^) is expressible as 



Gt''\t) 



I 

in '-^m+r-l,2m-l 

V 


\m"^ J 




^r — l,7n — 1 


[ — 1 

1 \m"^ J 


(1 - - - -)r"V 



or 



-,fc-n/^0,m+r-l 



^S''^ ((^) 



"-^r-l,m-l * 





/ 


Jk) \ 


1 


S 






\ 


^m,r J 


J_ 

m 




i yi — 1 



(1 



1 j \yn— 1 

?7i rn ' 1 
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Proof. From statement (i) in Theorem l2.1l it suffices to obtain the density of 
Smjr^kmlr ^ Pkm/r,n-km/r- Using the expression for f„i/r in (|3.3p . we note that 
a solution for the integral fjy^/^{tu){l — u)'^^^du, for g > 0, is obtained from 
Corollary 19. II as 

Q 1 iY^^ /• i— 1 \m 

\ ui 7^/1 'Vm-'l 

Q 1 J_\m-1 i j-l-q y 

V m m-'l 'V jn /I 



^m+r-l,2m-l 




Hence, absorbing the term t into the above expression with q = n ~ 

km/r by (|9.10p and taking the ratio of the resulting expression and (|3.3p gives 
'^m'/ri^y second result follows from simple algebra. 

We next show how to recover the result of Pitman (j40l . Corollary 4.11, p. 93) 
or (fsil ) in the case of m = 1 and r — 2, and show that the other rational cases 
of a may be expressed in terms of generalized hypergeometric functions. 



3.1. Hermite and generalized hypergeometric functions 



(i) When m — I and r = 2, we see that Theorem 13.11 recovers almost imme- 
diately the result of Pitman 

(0, 

Corollary 4.11, p. 93) or ((39l ) as follows. 



One has that 



1 (n,k) 
"1/2 



(t) is expressible as 





l+k k \ 
2 ' 2 \ 

— n J 














At 


:) 









(4i)-^-ie-*;7(-| 



it) 



(4t)- 



where U{a, h, c) is the confluent hypergeometric function of the second 
kind (see (j25l . p. 263)). The above ratio reduces to 



-k+l 



i"2 + 2{7 



k I 11 

h — , — 

2 2 ' 2 4t 



via an application of the recurrence relation (j47l p. 505) 
C/(a, 6, z) = z^-^U{l + a-b,2~b,z). 
A change of variable t yields the expression 2"-^'^ X'^^^ hk+i-2nW 



inside equation (110) in (|39[). where h^{X) is the Hermite function of index 
v (|25l . section 10.2), based on the following relationship, 



KiX) = 2-/^U 



2' 2' 2 
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(ii) Now, when m — 1 and r — S, G^j'^^'^\t) is expressible as 



G°;? I 27t 



2+k 1+fc 
' 3 ' 3 ■ 



G% I 27t 



1 _ 2 • 
3 ' 3 



where the G-functions at the numerator and at the denominator are re- 
spectively 



1^2(1 -n- 



fc. 1 2. 1 ^ 

3 ' 3 ' 3 ' 27t' 



1^2(1 -n+f + i; 



1.2 4. 1 

3' 3 ' 3 ' 27t 



r(n-|)r(i)r(|) 



3.^V3-V3. 3t^r(n-Mi)r(|)r(|) 

I 1-^ ^V-i^ " 3 ^ 3' 3' 3' 27t^ 



(3t^)2r(„_^)r(|)r(|) 



and 



27r 



0-^1 (' 3' 27t) _ (' 3' 27t) 

r(|) stini) 



wherein, for non- negative integers p and q, p < q or p = q + 1, \z\ < 1, and 
6, ^0,-l,-2,...j = l,...,g, 



p-Fq(ai 



ip;6i, ...,bq;z) 



[ai]j[a2]j 



[6l]^[62 



I ' 



(3.6) 



is called the generalized hypergeometric function (see (6, Chapter IV), (14£ 
p. 437)), which is available in Mathematica as 
HypergeometricPFQ [{al , . . . ,ap},{bl, . . . ,bq},z] . 

Remark 3.1. In general, when a = m/r with r > 2, both the numerator and 
the denominator of the ratios in Theorem \3.1\ are expressible in terms of sums 
of generalized hypergeometric functions of the form, 

2m — l-fm+r — 2 (ai, . . . , a2m-i; 6l, ■ . ■ , brn+r-2; ') 



L-Fr-2(ai, • ■ ■ , Qm-i; bi, . . . , br-2', ' 



respectively. 



4. Unconditional Gibbs 



Hereafter we assume that the mixing distribution can be represented as the 
density 

7„(i) =/l(<)/a(0 
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13 



Note that mixing the density ja over G^'''\t) leads to the foUowing class of 
operators 



Jo Jo 



(1- i-fc,i) 
(-n,l) 



dt 



which is a particular type of Fox iJ-integral transforms, see, for instance, (|2 
and W 



Theorem 4.1. Suppose that h(t) is a non-negative integrable function with 
respect to fa- Then, without loss of generality, set h(t)fa{t) — 1, and form 
the density "fait) — h{t)fa(t). Then the EPPF of PK{pa',"/a) random partition 
has Gihhs form 

k 

Pq,7„(?T-1, . . . ,nfc) = Vn,k a\nj-l 



where, for k — 1, . . . , n, 



Vr 



n.k 



a^-^V{k) 



Pka,n — ka 



s, 



a,{k—l)a 



l3{k-l)a+l,n-l-{k-l)a 



with Vi_i = 1. 

This result follows directly from Theorem 12. II 



4.I. G -transform for rational values 

When h{t) is set to be a G-function, one can use, for instance. Theorem 19. II in 
the appendix, to calculate many EPPF's as follows. 

Proposition 4.1. When a = m/r and h(t) is expressible as 



C X G:;% I at 
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where C is a constant and notation in Jff.jTp follows, then V„^k is 



cr(27r)'''(™-i) 



/ 



xG 



,(u+l)m+r— 1 
(a;+l)m+r-l,(ui+2)m-l 



{amy 



Tm^r, A(m, -di), . . . , A(to, -dx) 
(5m,r, A(m, -ci), . . . , A(m, -c„) 



where K^/r defined in i3.2\) . A{i,a), for any integer i, is defined in i9.19\) . 
and p and b* are defined for GJ^'^(crt) according to ^9.17^ . 

The result is just a specialization of Theorem 19. II given in the appendix. We 
now address two specific new examples. 



4-1.1. Example: Modified Bessel functions 

Proposition 4.2. Suppose we take 7i/2(^) oc Kjj{^/t)fi/2{t), or in other words, 
h{t) — Kri{-\/t)/ J Kjj{\/t)fi/2{t)dt, where Kr^{-) is the modified Bessel function 
of the third kind or Macdonald function (7, Sec. 7.2.1). Then, 



Proof. The result follows from Proposition 14.11 by substituting m — 1 and 
r — 1 and recognizing 





_1 


+fe k V V\ 






2 ' 2 ' 2 ' 2 \ 










~n ) 






1 »2 V\ 


^0,3 / 
^3,0 1 




2' 2' 2 \ 


16 





h{t) ^GxKniVt)^Gx Go ° 



2 ' 2 



by ([9T51) . where 

G-i = / K^{Vi)fi/2{t)dt 




G% I At 



dt 



due to Theorem 19.1 
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4- 1-2. Generalized Hypergeometric Functions as EPPF 



An EPPF in terms of one generalized hypergeometric function defined in (|3.6p (see 



Chapter IV)) results when Vn,k — /g ' [t)^{t)dt is representable as some 
constant multiplies either 



Gp.g+il-z), forp<g. 



or 



G^VtJ+i(-z), for \z\ < 1, 



(4.1) 
(4.2) 



due to JHHl). 

Proposition 4.3. Suppose that a = 1/r and h{t) — C x G^j\. I (r^) t 



where C is a constant, notation in {9.11^ follows, and none of (ci)'^ is identical 
to any of {dj)f. When 

(i) w > X + r, or 
(ii) X — w — r + 1 and w > r — 1, 



Vn,k is given by 



+ , - + I/, (d, + ; (cj + i^)5"; -1 



where v = 1 



CI 



i + k 



CK./rT^ IlZl r + ^) F ^) nr=i r(d. + ^) 



with Ki/j. defined in i3.2\) with m — \, provided that n ^ dj ^ 
Q ^ for i ^ 1, . . . ,w, j = 1, . . . , X, £ ^ 1, . . . ,r. 



k+e-i 



Proof. It follows from Proposition 14.11 with m = 1 that Vn,k is equal to 
CKi/^r^ multiplies 



G 



l,x-\-r 
x~\-r,w-\-l 



f^l,x+r I -, 
^x+r,w + l ^ 



( i+k \r~l k ( 

V r 'I ' ill 



(1 - (i±fe + ^));-\ 1 - (I + ^) , (1 - (d„ + ^))f 
0, (1 - (c, + v))^ 



followed from (|9.10p . The last G- function takes either the form of (|4.1[) when 
ui > a:: + r, or the form of (14. 2|) when x = w — r -\- 1 and w > r — 1, and, hence, 
the result follows from ()9.16|) . 
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5. First distribution theory example: Sa,9 given Xa,e 

We now come to our first result which does not rely on special properties of 
G- or iJ-functions and importantly applies to all values of a. Let Sa and Sa,9 
denote independent random variables having laws described previously. Then, 
define the random variables 

Xo.,e = 1^ (5.1) 

whose laws have been recently studied in James f|18|). They represent a natural 
generalization of the random variable 

where Xa = X^fi and S'^ is independent of Sa and has the same distribution. 
Remarkably although Sa does not have a siinple density, except for a — 1/2, 
Lamperti l|24r ) (see also Chaumont and Yor |4|, exercise 4.2.1)) shows that the 
density of Xa is 



sin(7ra) y" ^ 

y2a _^ 2?/"cos(7rQ;) + 1' 



fx. {V) = : ' oa . , for y > 0, (5.2) 



with cdf 

1 / v° 

Fx. (y) = 1 cot^i cot(7ra) + ^ 



sin(7ra) 

Note furthermore, when a = 1/2, 



d ^1/2 , ^ d ^9+1/2 

- — and Xi/2,e = 

^1/2 ""1/2 



^ d "i/2 , ^ d "«+i/2 .J. „x 

= 77^ and Xi/2,e = ^ (5.3) 



where G'j^yj — ^"1/2 and all the gamma random variables are independent. 
See (firl . section 4.2) for more on the variables (|5.3p . 

Here we investigate the mixing distribution 7 corresponding to the law of 
Sa,e given Xafi = 1- That is, the random variable with density 

la{t)^Ca,et'-'fa{t)fa{t) (5.4) 

where 



So from (E 



Ca,e^l/fx.A^). 
27r(l + cos(7ra)) 



Gq — 



sin(7ra) 



Now define, for 6* > 0, 



■ 2x" cos(a7r) + 1] 2° 
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We now give a brief description of the density of Xa,e for > and general a 
which is due to James (l8;, Theorem 3.1). When 9 = 1, 

f I \ \ ( \ ^ sin(7rFx„(y)) 



TT \ry2a _|_ 2ya cos(a7r) + 1] 2" 

Hence, the normaUzing constant in this case satisfies 

[2(l + cos(a7r))]^ 

In general, for 6* > 0, 

fx^,e{y)= r iy-xf-'A'gix)dx (5.5) 
Jo 

where, suppressing dependence on Fx^, A' denotes the derivative of A. Fur- 
thermore, define 

^(„,fc) ^ [sin(7ra)]'''^" (x" cos(a7r) + l)^x"(''-J) 
J TT [2;2a _^ 2a;« cos(a7r) + 1]'' ' 

Before we proceed with the description of the EPPF, we provide a representation 



Lemma 5.1. For < a < 1, and fc = 1, 2, , 

k 



Proof. Since k = 1,2,..., we first apply the multiple angle formula to 

sin(7rfcQ!Fx„ (x)). 

The result is concluded by noting the following identities which are given in 
James (fisl . Proposition 2.1), 

sH^aFxJx))- 



+ 2x« cos(a7r) + 1]^^^ 
and 

x" C0s(Q;7r) + 1 

cos[nal:<x^(x)) = tj-. 

[a;2" + 2x°' cos(a7r) + 1] ' 

Then we have the general description of the EPPF. 



Proposition 5.1. Suppose ja is specified as |5.^| j for 9 > —a. Then, for n 
2,3..., and k = I, . . . ,n, 

a'^-'T{9/a + k) .n+e-2 , 



Vn,k = C^.e ' / (1 - xY+'-'Ae+ko.{x\FxJdx 

1 (fc* + n — 1) 
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(i) In particular, when 9 = 0, we obtain 



a*-ir(/c) ^ fk 
r{n- 1) 



{n,k) 



(5.7) 



whe 



■' sin(7ra) Jq ■' 



Note that df'^\x) > /or < a; < 1. 
(ii) When 9 = and n > 1, 



Vn.l = 



2(1 + cos(a7r)) (1 - a;) 



T{n - 1) Jq [x^" + 2x" cos(Q;7r) + 1] 
Proof. It is easy enough to work directly with the expression 



dx. 



(5.8) 



r(n - ka)fa{t) 
Mixing relative to 7a (0 defined in (|5.4p . one sees that 



fa{tu){\~u)^ ^du 



,i-{e+ka), u..,, r{9/a + k + l) 



^Utu)Ut)dt = r(g + fc^ + i) /^-.«+--(")- 



Now noting the form of the density of Xa.e+ka from (|5.5p and integrating with 
respect to u lead to the evaluation of the integral 



{u-xy+'^'^'ii-u) 



i—ka—l 



du = (1 — x) 



n+e-i r(6' + ka)T{n — ka) 



T{9 + n) 

Ignoring constants for a moment this leads to an integral of the form 



1 - x) 



Now, since 9 > —a, it follows that n + 9 > 1 when ti = 2, 3, . . . . This allows us 
to use integration by parts to get 



{l-xr+''-'A',+,^ix)dx = in + 9-l) I il^xr^''-'Ae+ko.{x\FxJdx 



which yields the general expression for 9 > —a. Statements (i) and (ii) then 
follow from an application of Lemma 15. li and also the use of (|5.2p . 

When a = m/r, we may easily express Vn,k in terms of G- functions, which 
we leave to the reader. We focus on two interesting cases. 
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5.1. 2F1 EPPF 

When a = 1/2, we have that for 9 > -1/2, 



Ci/2,e — 



r(0 + i)" 



One can then show that 



k 

2 ' ' 



k , 1 
2^2 



0,-61- n + i 

Now applying (|9.16p or identity (2.9.15) in (20), the last G-function reduces to 



r(0 + f + i)r(0 + f + i) 



T{0 + n+^) 



k 1 

2 + 2^' 



n+-;-l 



where 



2F1 {a,b;c;x) = ^ 



r(c) 



[c], ~ r(6)r(c-6) 



(1 — a;r)" 



is the Gauss hypergeometric function (see (jSJ)), which is a special case of 
Note this result can also be checked by using the explicit density of 



5.2. 3F2 EPPF 



When a — 1/3, that is, m = 1 and r = 3, one gets 



^39 



Gli 1 1 



[r(i/3)r(2/3)]2 3,3 

Furthermore, the G-function reduces to 

r(g + | + i)ntir(g + ! + i) 
r(i)r(|)r(n + + i)r(2(0 + ^)) 



k _ 2 _k _ 1 _k 
3 3' 3 3' 3 



3-P2 



fc+2 „ _ fc 
3 ' 3 ' 



fc+3 
3 



^+ 1 2(^+fe±3) 



6. A recursive method for calculating EPPF's and Fox ff-transforms 



Definition 3 or equation (8) of Gncdin and Pitman (|l2f ) establishes the following 
backward recursion for all Vn.k, n — 1, 2, . . . ; fc = 1, 2, . . . , n. 



Vn,k = {n- ka)Vn+iM + K+i,fc. 



fc+i 



(6.1) 



with Vi^i — 1. 
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One can turn (|6.1|) around to obtain the foUowing forward recursion, 

Vn+l,k+l = VnM - [n- ka)Vn+i,k- (6.2) 



The key point about the recursion ()6.2p is that it enables computations of all 
Vn k for n = 1, 2, . . . ; A; = 1, 2, . . . , n, provided that 



Vn.l = 



1 



Tin) 



E 



Sa 



(3l,n-l 



for n — 1,2,..., are known. For clarity, we demonstrate this point by considering 
the case of n = 3. When Vi^i, ¥2,1 and ¥3^1 are given, one can obtain all Vn,k for 
n = 1, 2, 3, A; = 1, . . . , 71 as follows. 

(i) compute ¥2,2 — Vij — (1 — a)V2,i; 

(ii) compute ¥3,2 — ¥2,1 — {2 — ce)¥2,2', and 

(iii) compute ^3^3 = ¥2,2 - (2 - 2a)¥3^2- 

So quite fortunately we can focus on the relatively simpler task of calculating 
¥n.i- The simplicity occurs because the quantities only depend on the distri- 
bution of the independent pairs (S'q, /3i_„_i), or equivalently, {Sa,a, Pa,n-a)- In 
particular, when Sa does not depend on k or n, we show some interesting appli- 
cations of the recursion in section [HI One may also see how this applies to the 
explicit expressions in (|5.8|) and (|5.7|) . 



6.1. An all purpose solution? 

It is widely believed that beyond the infinite series representation of fa, there 
are no general explicit representations of fa ■ In fact this is false as one may use 
the representation of Kanter (fiol). which we now describe. Setting 



Kaiu) 



sin(7rau) 


1 

l-Q 


sin((l — Q!)7ru) 


sin(7rit) 




sin(7rau) 



it follows from Kanter (see also Devroye (5) and Zolotarev ((sol)) that 

3-""^^°(")K„(M)du. (6.3) 



-1/(1-") 



1 — a 



That is, 



Sa - (K„(t/)/Gi) 



(l~a)/a 



(6.4) 



where U is a uniform random variable on [0,1], independent of a standard 
exponential variable Gi. So, using (|6.4p . one can obtain calculations for V„,i 
and hence 14,, fc by simulating U,Gi and a variable. In addition, one can 

certainly obtain new integral representations for ¥n,i and ¥n,k- 

However, while this is certainly true in principle, and useful in some cases, 
there is still quite a bit lacking in this representation from an analytic point 
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of view. For instance, it is not obvious how to use (|6.3|1 to obtain the Laplace 
transform of Sa or to obtain the nice form of the density of Xa- In addition, 
one would require that h(t) has a manageable form, which is not the case for 
instance in section [8] Recall that even the known tractable case of a = 1/2 does 
not always immediately lead to nice expressions for We will see in the next 
sections that one can do quite nicely without resorting to (j6.3p . 



7. Classes of EPPF's via Probability Transforms 

The first sections focused on the fact that we could use the fact that Mejier 
G-transforms can now be easily computed to calculate Vn,k and related quan- 
tities. We also noted that while in theory the Vn,k can be calculated based on 
7?-transforms, the practical tools to do this are not yet available as not all H- 
functions can be computed easily without resorting to numerical integrations. 
We now demonstrate through using probability distribution theory how to cal- 
culate Vn,i and hence by the recursion all {Vn,k)- This in turn will provide new 
methods to calculate many i?-transforms as well as G-transforms. We will focus 
on three general classes, but certainly more can be constructed. 



7.1. Lamperti Class 



Here, we introduce an entire class of EPPF's based on Xa and Xa.ka- First, for 
a positive integrable function g, let 



1/La 



sin(7ra) 



g{y)y°' ^dy 



Ti" Jo + 2?/" cos(7rQ!) + 1 



-%(^a)]. 



Proposition 7.1. Suppose that for each < a < 1, h{t) = LaE[g{S'^/t)]. That 
is, 

7aW=iaIE[5(5;A)]/a(i)- 

Then, PK(pa,Ja) has the following properties. 
1 



ft) Vn 



Tin) 

fti) For k — 1,2, ... ,n, 



LaE[giXaPi,n-i)]- 



„fe-i 



Vnk — 



m 



r{n) 



LaE[g{Xa,kaf3ka,n-ka)]- 



ftii) Equivalently, 



r(?i) 



5 |/3.,„-fc^X(^)z?, 



where {xi^^) are iid random variables equal in distribution to Xa and, in- 
dependently, {Di, . . . Dk) is a k-variate Dirichlet{\, . . . , 1) random vector. 
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Proof. Statements (i) and (ii) follow from Theorem 12.11 and the definition, 
Statement (iii) amounts to showing that 

k 

-^ajkaPkotjU — ka — Pk,n—k ^ ^ -^a^-^j' 



In order to do this, we refer back to the recent work of James dlilB). First, 
write ^ 

Pka,n—ka — f^k,n — kPka,k{l — a)- 

Then using the notation and the result in (18, Proposition 3.1(ii)), X^^ka = 
Mka{Fxa)- That is, X^ka is a Dirichlet mean random variable of order ka 
indexed by Fx^ ■ Let denote a Bernoulli random variable with success prob- 



ability a, independent of Xa- Now, applying Theorem 2.1 in James (llTI ). it 
follows that 

/3fca,fc(l-a)^fca(-Fx„) = Mk{Fx^YS} 

where Mk{Fx„Yj) denotes a Dirichlet mean of order k indexed by the cdf of 
XqFq,, Fx^Ya ■ Again, these specific random variables are found in |1^ ). Now an 
application of Proposition 9 in Hjort and Ongaro (fl^ (see also Proposition 4.5 
in (0)) yields 

k 

MkiFx^Yj = D,m[^\Fx^yJ 
where m\^\Fx„y^.) = Mi{Fx^y^) are independent. Now from Proposition 



3.1 (iii) in James (|l8r). one has 

X^ = Mi{Fx^yJ, 

which completes the result. 

In view of this result, we proceed to give an expression for the density of 

-^a,kcxPka,n—ka- Define 

/•1/w 

/„„^ — /"^ 1^ / ^^ Q(n,fc) , 



iPp^' {w) = {n-l) (1 - wx)"-^^)"'^' {x)dx 

Jo 

where ??^"'''^(a;) is given in (|5.6p . 

Lemma 7.1. For n = 2, 3, . . . and k — 1,2, ... , the density of Xa.kaPka.n-ka 
in w is given by 
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which, from Lemma \5.1\ can also he expressed as 

3 



3=0 



Proof. Using standard arguments one can write the density of Xa,kaPka.n-ka 
as 

rr^y°l ^ - y^lyT-'''-'v-'Vx^.,Ay)dy. (7.1) 

1 {ka)i [n — ka) 

But it is not difficult to show that for any 6 > —a, 

Hence making the change r = 1/y, the expression in (|7.ip is a constant multiphes 



Now arguing as in the proof of Lemma l5.ll based on the definition of the density 
of Xa,ka provided by (j5.5p . we will arrive at the integral 

vj-^ f\l-u)"-'''"\u-tw)''"-^du 



where t refers to the argument in A^^(<). So, it follows that the density in (|7.1 
can be expressed as 

- f"{i-wtr-^A',^{t)dt, 



w 

and the result follows. 

This leads to another description of the quantities in Proposition 17. II 

Proposition 7.2. Suppose that for each < a < 1, 

Jait) = L^n9{S'Jt)]f^it). 
Then, the EPPF for PK{pa,^a) is such that, 

where 



.(„,fc) _ [sin(7ra)]" ' p ^ ( l3i,n-iW (rcos(7ra) + 1)' r^^'^)-^ 



„l/a 



[r^ + 2r COS {air) + 1]'^ 



dr. 



Proof. This follows from Lemma mi and Proposition l7.1l and applying a change 
of variable. 
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7.1.1. Example: Mittag Leffier and generalized Mittag Leffier functions 



This example is also influenced by some results in (|l8r ) and concerns the Mittag 
LefHer function and some of its generalizations. Recall that the Mittag Leffler 
function may be defined as 



1=0 



{l + la) 



Efe 



Furthermore, there is an integral representation 

sin(7rQ:) 



E„,i(-A) 



_^ 2?/" cos(7ra) + 1 



dy 



which allows one to compute the Mittag Leffier function. Now, define 



M)' [fc + i], 

l\ T{\ + ka + la)' 



Theorem 7.1 in James (18) shows that 



(7.3) 



(7.4) 



In this section set L^.x — 1/Eq^i(— A). 



Proposition 7.3. Let, for < a < 1, (7q,a : A > 0) denote the family of 
densities each defined as 

7a,A(i) =ia,Ae"^/*°/a(t). 

Then, the EPPF's of the PK(pa,^a.\) family satisfy, for each fixed (a. A), 

ft) VnA = f^La,AlE[E„,i(-(/3i,„-i)"A)]. 
r(n) 

(ii) For k — 1,2, ... ,n, 



Vn,k 



)"A)]. 



(Hi) Equivalently, 



Vnk — 



m 



r{n-k) Jo 



Y[E^A-pTb''\)dpi 



I' 1=1 



b'^-'il-b)' 



where Sk = {{pi, . • . ,Pk) ■ < J2i=iPi < !}• 
(iv) Vn,k is described by Proposition\7.S\ with 



jn,k) _ [sin(7rQ;)]'' ^ 



i/°/3i.„_i/r^/°j ('^'^os(7ra) + I)-' ^ 
[r^ + 2r cos [air) + 1]'^ 



dr. 
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Proof. First note that the relevant g function is 

5A(x)=e-^''°^ 

Hence statement (i) follows from statement (i) in Proposition 17.11 and the rep- 
resentation of the Mittag LefFler function. Statement (ii) follows from (|7.4p and 
Proposition 17. II Statement (iii) follows from statement (iii) in Proposition 17. II 
and statement (iv) follows from Proposition 17.21 S tatement (iii) is also a special 
variation of statement (v) in Theorem 7.1 in (jlST ). 



7.2. Beta- Gamma Classes 



The previous section identifies a large class of EPPF's that can be computed 
by using the random variables Xa and Xa,ka- Compared with common random 
variables that one encounters in standard probability textbooks, these random 
variables are rather exotic. In this section we show, also using results that appear 
in (|l8l). that we may define a large, and clearly important, class where the 
EPPF's are expressible in terms of expectations depending only on beta and 
gamma random variables. 

Recall the identity in (|1.5p . which again may be found in Pitman (fiol section 
4.2) and Perman, Pitman and Yor |35l ). Statement (iii) in Proposition 3.2 in 
James (fisl ) uses the result to establish the identity 



Gg 



+1 



(7.5) 



for 9 > —a. James (|18l . section 3.0.1) also explains how this amounts to a 
rephrasing of an otherwise equivalent result in Lemma 6 in Bertoin and Yor (0). 
Now let 

l/I]„,, = E[5(Ge/.+i/3eV,i-a)]- 
Then (|7.5p suggests the following class. 

Proposition 7.4. Suppose that for each < a < 1, and 9 > —a, 



Ga 
0+a 



That is, 



hit) 



T{9 + 1) 



T{6/a + l) 

Then, PK{pa,"fa) has the following properties 
(i) For n= 1,2,..., 

T{9+1) 



Ga 



r(n 



^aMMgiGe/a+lf^e+a.n-a)]- 
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(ii) For k = 1, 2, . . . , n, 



a 



fe-i 



r{e + i)T{0/a + k), 



Vn,k - r(6i /a + l)r(n + 9) ^"-^^If ('^«/"+'=^e+a,n-a)]- 
Proof. It suffices to examine the quantity 

ka,n — ka }\ 

for each k, g and fixed y. Now noting the form of the density of Sa.e and that 
of a beta random variable, it follows that the above expectation can be written 
as 

r(n)r(g/a + fc) 

r(n + 9)r{k) ^S^y afi+kaPe+ka,n-ka)\- 

Now replacing y with Gg_|_^ shows that 

_ a'--ir(g + i)r(0/a + fc) ^ 

•^".fe — r(0 /a + l)r(?7 + 0) "'^ L9l^6l+Q'^a,0+fcaPe+fcQ,n-fcQjJ- 

But, using the calculus of beta and gamma random variables, it follows that 
since 6* + a < 6* + /ca, 

Ge+a r, Ge+ka 

Pe+a,(k-l)a- 



iJa,e+ka Oafi+ka 

Additionally, 

The result is concluded by applying the identity in (j7.5p . 

Note that one can also work with the other equality in (|7.5p . We do not 
discuss that here. 



7.2.1. Example: Hermite Type 
Note that 

1 r°° 

E[e-G?+»/*°] = —1— / e-(f )°a;^+«-ie--da: 
T(e + a) Jo 

which looks like some sort of generalized Hermite function. We can use this fact 
along with Proposition l7.4l to calculate an EPPF, based on the mixing density, 
which we deliberately write in the less obvious form, 

7o.A W = ^-'^l^^l^ f r e-^(f )°x«+"-^e-dx (7.6) 
r(6' + a) Jo 

where 

T{e + a)T{l - a) Jq 
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which is the same as 
Here again 

g{x) = e"^, 

which leads to the next resuh. 

Proposition 7.5. Let ja,\ be specified by ( [7.6'[ j for all < a < 1 and 9 > —a. 

Then, the EPPF of PK{pa,ja,\) is specified by 

= $W (A) /V + A^.")-(«/"+^)^.''+"-i(l - ur-"-'du 

where 







(fe) ^ a^-ir(g + l)r(g/Q + fc)S]^.e(A) 
"^^^ ^ T{e/a + l)T{9 + a)T{n-a) ' 

7.3. Comment on Composition Classes 

Recall that if one evaluates a stable subordinator at a positive random time, 
say T, it is equivalent in distribution to 

having Laplace transform 

where e"'^^'^^ is the Laplace transform of T evaluated at A. We see that the two 
classes in the previous sections are special cases of this where one chooses T in 
such a way where we can work more easily with the stable random variable. It 
is clear that many models may be derived using this general idea. That is, based 
on 

/i(t) = E[5(tri/")]. 

Furthermore, since g is arbitrary one can create quite a few variations of this. 

8. EPPF's via ranked functionals of self-similar Markovian 
excursions 

We close this paper by looking at a large class of mixing densities that probably 
would not immediately come to mind. This is partly because, as we shall show, 
the apparently complex form of them. The mixing distributions are based on 



the results of the interesting paper of Pitman and Yor (j44f ). In fact, in the 
general case the densities do not have known closed forms, which means that 
h{t) does not have an explicit form and hence brute force simulation methods 
or the calculus of H- or G- functions cannot be used directly. However, as we 
shall show, these are members of the Beta-Gamma class, and we will indeed be 



able to get explicit results by exploiting Corollary 7 in Pitman and Yor (|44[ ). 
For clarity, we start off with a simpler yet still challenging case. 
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8.1. Example: Hyperbolic Tangent and Kolmogorov's Formula 

First note the famous random variable 

M"^^ := max iS^n 

0<u<l 

where Sjj'' denotes standard Brownian bridge on [0,1]. One has the Kolmogorov's 
formula, 



1^-21' 



mil' < ^) - E (-i: 

Furthermore, one has the remarkable formula 

n\Bi\M^' <y)=tanh{y), (8.1) 

where Bi is a standard Gaussian random variable independent of B^'^ . See Pit- 
man (|40l . p. 203) for this description. We will try to be quite transparent at this 
point. We want to use Proposition 17.41 to obtain nice EPPF's generated by a 
density related to A/|"'. The key to exploiting these results is that 

Looking at Proposition l7.4l we want this to apply for all a, we solve the equation 



- + 1 = 1/2 
a 

which leads to 

= -a/2, 

which is less that but greater than —a, hence in the acceptable range of 
Proposition l7.4l This leads us to construct the following mixing density for each 
T > 0, 

7a,a/2.rW = /a,-a/2 WP(2(M^'0' < rr/G^/2)Sa,a/2 M (8.2) 

where by Proposition 17.41 and (|8.ip . 

l/So,a/2(r) = Pm\M^^ < V^P^^ilj = E[tanh(V^/3;;,/J_J]. 
Now define probabilities 

PnAr\a) = mCk^iMM^'f < K%,u-o.-')- (8-3) 
This leads to the following result. 

Proposition 8.1. Let j or < a < 1, and each r > 0, PK{paT"ia.a/2,T) denote 
the Poisson Kingman distribution with mixing density specified in \8.2]) . In addi- 
tion, letpn,k{T\a) denote the probabilities defined in i8.S\) . Then, PK{pa,ja,a/2,r) 
has the following properties. 
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(i) For each n, 



r(i - a/2) n^^MV^P^jll-J 

VnA — 



r(n-«/2)E[tanh(V?/3:/tLj]' 
(ii) For k — 1,2, ... ,n, 

a'=-ir(l-a/2)r(fc-l/2) p„,fc(r|a) 



Vn 



k 



r(l/2)r(n - a/2) E[tanh(V^/3;;(J_ J] ' 

(Hi) For T > 0, applying \6.^) yields the following recursion, 

(n — a/2) 

Pn+l.k+l{T\a) = -J^[Pn,k{T) -_p„+l,fc(T)] + P„+l,fc(-r). 

(fc — l/2)a 

Hence, these probabilities are completely determined by 
E[tanh(V7/3;;2/tj]=Pn,i(r). 

Proof. For clarity, it follows from the form of the density in (|8.2p that Propo- 
sition [7]3] applies with 

gix) ^ P{2{M]")^ < t/x) 
and d = —a/2. Furthermore, 

E[5(Gfc_i/2/32/2,n-a)] = VnAA^) 

which specializes in the case of fc = 1 to 
due to 



8.2. General result 

Next we first sketch out some of the notation and definitions given in Pitman 
and Yor (44). We ask the interested to consult that work for more details. 
Note instead of 7 and a, used in that work we use k and S. The first change 
is obviously to avoid confusion with the mixing density. The second as in the 
previous result is that S will be associated with a stable subordinator, and other 
relevant quantities, not directly related to our use of a. We now paraphrase 
Pitman and Yor (|44 p. 366-367) which describes the general random variables 
we are interested in. 

In this section B := {Bt,t > 0) denotes a real or vector valued /3-self similar 
strong Markov process, with starting state which is a recurrent point for B. 
Let (et, < t < Ve) denote a generic excursion path, where Ve is the lifetime or 
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length of e. Let be a non- negative measurable functional of excursions e, and 
let K> 0. Call F a n-homogeneous functional of excursions of B if 

F{euQ< t<V,)^ V,^F{V,'f'e^v.,Q< u < 1). (8.4) 

In particular, we (Pitman and Yor) have in mind the following functionals F: 
length, maximum height, maximum absolute height, area, maximum local time. 
Now we restate their Theorem 3. 

Theorem 8.1 (Pitman and Yor [a^). Let F be a n-homogeneous functional of 
excursions of B , let F^ :— F{B°^) for a standard excursion B^^, and suppose 
that E[(i<"*)*/''] < oo, for < 5 < 1. Then the strictly positive values of Fie) as 
e ranges over the countable collection of excursions of B^'^ can be arranged as a 
sequence 

F^' > F^' > ■■■ > 0. 
Let Ts be a random variable, independent of B^'^ , with the gamma{5, 1) density 

F{Ts £ dt)/dt = T{5y^t^-^e-^ for t > 0. (8.5) 

Fix A > 0. Then, the joint distribution of the sequence {Fj"-',j — 1,2,...) is 
uniquely determined by the equality in distribution 

(M(A-^r?Ff ), j = 1,2, . . .) ^ {T*,j = 1,2, . . .) where T* := ^e./eo (8.6) 

1=1 

for independent standard exponential variables eo,ei,£2, • ■ ■, o,nd fj, is the func- 
tion determined as follows by S,X,k and the distribution of F^,: 

M^):-y^ j:^^t-'-^e-^*F{F^>xt-'^)dt. (8.7) 

As noted by Pitman and Yor (44) , the above result gives a characterization of 
the random variables Fj'^' but it does not provide an explicit expression for their 
distribution, as for instance in the special case of Af}"'. Nonetheless, letting 
denote the law of these F^"-', depending on S as described above, we can construct 
a rich class of EPPF's from these random variables. Specifically, Theorem 3 and 
Corollary 7 in Pitman and Yor ('S) along with Proposition 17.41 suggests the 
following construction of mixing distributions. 

First we solve the equation 

a 

which leads to 

-a < 9 = {S ~ l)a < 0. 
Now construct mixing densities as follows. For each a, j and S, 

7a,,,5(t) = (S,) X {X-'^F^^ > wit/GscT'-) /„,(5-i)a(t) (8.8) 
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where 



and /Lt defined in (18.71) is a quite general function, some examples are given in 
Pitman and Yor (HJ). Now define 



-CXK 



(8.9) 



The next result follows from the above discussion combined with Proposition l7.4l 
and a bit of algebra. 

Theorem 8.2. Let, for Q<a<\, Q<5<\, j — 1,2,..., PK{pa,^a,j,5) 
denote the Poisson Kingman distribution with mixing density specified in 18.8\) . 
In addition, let pn^k{w\S, a, j) denote the probabilities defined in i8.9\) . Then, 
PK(pa,"fa,j,s) has the following properties. 

(i) For each n. 



(a) For k — 1,2, ... ,n. 



E < 

r(l- (l-(5)a) 










T{n-{l-5)a) 

E< 






'1 


+ 



a'=-ir(l - (1 - 6)a)V{k - (1 - 5)) 


Pn,k{'w\5,a,j) 




T{5)V{n^{\-5)a) ( 

E <^ 




] 







(Hi) For A > 0, applying h6.2fl yields the following recursion, 

{n-{l-5)a) 



Pn+i,fe+i(w|Q;, (5, j) 



{k + 5 - l)a 

+Pn+i,k{w\a,5,j). 



'Pn,k{w\0L, 6, j) - p„+i^fe(w7|a, S, j)] 



Hence, these probabilities are completely determined by 



i^i^pas,:-j 

1 + 



Pn,iiw\a,S,j). 
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Proof. The result again follows from Proposition 17.41 Now from (|8.8p . fixing 
j, a, S and w, 

g{x) ^FsiX^^F^' > wx-"). 
Furthermore, 9 ~ {S — l)a. Hence from Proposition 17. 4| 



E[g{Gk+S-lPsa,n-a)] = Pn,fc(w|5, j) 

and when fc = 1, one uses Corollary 7 or Theorem 3 in Pitman and Yor lii) to 
obtain 



1 + 



1 3 



8.2.1. Bessel Bridges 

As an interesting special case, we now follow the example given in Pitman and 
Yor (0, p. 375). Let 



h-s{x) = 



denote the ratio of two modified Bessel functions, which follows as a special case 

of ([87]). Then setting k ^ 1/2, A ^ 1/2, Fj" ^ M]'^ where under Ps are the 
ranked heights of excursion of a standard Bessel bridge of dimension 2 — 26. 
Hence, in this case, 

g{x) =Fs{V2M^' > w2;-i/2), 
and we will use, from Pitman and Yor (|4^ . p. 375), 



Thus Theorem 18.21 specializes as follows. 

Proposition 8.2. Consider the setting in Theorem \8.Si with X — 1/2 and n = 

1/2, such that for each j, F^'^ = Mj^ are the ranked heights of excursion of a 
standard Bessel bridge of dimension 2 — 26. Then, 



and the Vn,k «re obtained from Vn,i via the recursion formula based on obvious 
adjustments in Theorem \8.'A This result generalizes Proposition \8 . 1\ 



r(i 


-(1 




r 




'} 


T{n 


-(1 






'l-h.siwl3jiD 


'} 
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9. Appendix: Fox H- and Meijer G-functions 



Definition 9.1 Function). For integers m,n,p,q such that < m < q,0 < 
n < p, for ai,bj e C with <C, the set of complex numbers, and for ai,l3j G 
K+ = (0, oo), i = 1, 2, ... ,p; J = 1 , 2, . . . the H -function is defined via 
a Mellin-Barnes type integral (see fSa: IH; \3sl: \2^) } on the complex plane in the 
form 



Tj7n,n 



1 

27ri 











(ai,ai), . 




z 




Tjm.n 

= ^P-q 


Z 












(&i,/3i),. 





llr{b,+f3,s)Y[Til^a,-a,s) 

P Q 

Yl r(a, + a,s) Yl Til-bj-P,s) 

i—n-\-l J— m+1 



-z-'ds. (9.1) 



Here 



z = exp [— s(log |z| + i argz)] , z 7^ 0, 



1, 



where log|z| represents the natural logarithm of \z\ and argz is not necessarily 
the principal value. An empty product in i9.1]) . if it occurs, is taken to one, and 
the poles 



b.ji^ — ^ , j = l,...,m;/ = 0,l,2,. 

Pi 

of the gamma functions T{bj + [3js) and the poles 

—Oi — k 



aik = 



i — 1^ ... ,n]k — 0,1,2, ... . 



(9.2) 



(9.3) 



of the gamma functions T{ai + ons) do not coincide: 

ai{bj + I) ^ f3j{ai - k - 1), i ^ 1, ... ,n;j ^ 1, ... ,m;k,l = 0,1,2, 

£ in 119.1]} is the infinite contour which separates all the poles bji in \9.'2\) to the 
left and all the poles in i9.S\) to the right of and has one of the following 
forms: 

(i) £, = £_oo is a left loop situated in a horizontal strip starting at the point 
—00 + iipi and terminating at the point —00 + iip2 with —00 < tpi < ip2 < 
+00; 

(ii) £, = £+00 is a right loop situated in a horizontal strip starting at the point 
+00 + iifi and terminating at the point +00 + iip2 with —00 <ifi<ip2< 
+00; 

(Hi) £, = £i-yoo is a contour starting at the point 7 — ioo and terminating at the 
point 7 + zoo, where 7 e M = (—00, 00). 
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Refer to Kilbas and Saigo (|20l ) for more discussion about _ff-functions, and 
in particular, Theorem 1.1 for the situation in which H™^"{z) defined by (|9.ip 
makes sense. 



Definition 9.2 (Meijer G-Function). In general the Meijer G-function l( 291:13 d: 
31) is defined by the following Mellin-Barnes type integral on the complex plane: 



1 

27ri 



P,9 



ai, 
bi, 



,bq 



l[Tib,+s)Y[ril-a,-s) 

P Q 

n naj+s) n r(i-6,-s) 

j— 71-^1 j—m~\~l 



-z-'ds, (9.4) 



where the contour of integration £ is set up to lie between the poles of T(ai + 
s) and the poles of T(bj + s). The G-function is defined under the following 
hypothesis 

• 0<m<(7, 0<n<p and p < q — I; 
. z ^ 0; 

• no couple ofbj, j = 1,2, ... ,m differs by an integer or a zero; 

• the parameters G C and bj G C are so that no pole of r(6j + s),j = 
1,2, ... ,m coincide with any pole of T{ai + s),i = 1,2, ... ,n; 

• Oi ~ bj 7^ 1,2,3,... for i = 1,2, ... ,n and j = 1,2, . . . , m; and 

• if p ^ q, then the definition makes sense only for \z\ < 1. 

Refer to (@, Sec. 5) for a more thorough discussion of the G-function. Both H- 
functions and G-functions are very general functions whose special cases cover 
most of the mathematical functions such as the trigonometric functions, Bessel 
functions and generalized hypergeometric functions. Nonetheless, G-functions, 
but not //-functions, are implementable in Mathematica as Mei jerG [{al , . . . , an} , {a(n+l) . 
{b(m+l) , . . . ,bq},z]. Comparison between definitions ()9.1|) and ()9.2|) reveals 
that any G-function is an //-function, but not vice versa; when ai = Pj — 1 for 
i = 1,. . . ,p and j = 1, . . . ,(?, 



. ,ap},{bl, 





(a.,i)r 






z 




— (~im,n I 

— Pi9 1 











(9.5) 



This implies that those //-functions which reduce to G-functions can be evalu- 
ated in Mathematica. Here are some results, available from Kilbas and Saigo ([2~ 
Section 2), that follow directly from the definition of the //-function in 



Property 9.1. The //-function is symmetric in the set of pairs 
(ai, ai), . . . , (a„, q;„); in (a„+i, a^+i), . . . , (cp, a^); in {bx,(i\), . . . , (6m, /3m) and 

in (fem+l,/3m-|-l), ■ • ■ , {bq,Pq). 
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Property 9.2. If one of the pairs {ai,ai),i = 1,2, ... ,n, equals one of the pairs 
{bj,pj),j = m + 1, . . . ,q, and n> l,q> m, then, for instance, 



(ai,Q!j)i 

{ai,ai) 



TTm,n — l 



(9.6) 



or, if one of the pairs {ai,ai),i = n + 1. . . . ,p, equals one of the pairs {bj,Pj),j = 
1,2, ... ,jn and rn > l,p > n, then, for instance, 



Trm—l,n 



Property 9.3. There holds the relation 

(ai,ai)i 



p-i 



Tjm.n 



ibj,Pj)l\ 



Tjn.m 



(1 - aj,aj){ 



Property 9.4. For y > 0, there holds the relation 













{ai,kai)\ 


TTm,n 

p,g 


z 














{b,,P,n. 







Property 9.5. For cr G C, there holds the relation 

{tti + aai,ai){ 

(6,+a/3,-,/3,)? 
Property 9.6. For a,b G C, there holds the relations 



^ ^p,g 





(ai,ai)? 






z 




= ^P,g 


z 




{bj,Pj)l_ 





^P,g 



(a,0), {ai,ai)l 



{bj,Pj)l 

when Re(l — a) > and n> 1; 

jjm,n ^ 

[ {b,Q),(b„P,)l 
when Re(6) > and rn > 1; 

K,ai)i~\(a,0) 

Tjm.n 

^P,g 

{bj,m 



V{l-a)H^jJl 



m,n— 1 



r(6)//- 



m — l,n 
9-1 







z 






ibj,/3j)l_ 






z 






{bj,Pj)l 



l>) p-^'" 



(bj,P,)l 



(9.7) 



(9.8) 



(9.9) 



(9.10) 



(9.11) 



(9.12) 



(9.13) 
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when Re(a) > and p > n; and 

{b,,p,)r\{b,o)_ 

') > and q > m. 



Tjm.n 



r(i -6) 



9-1 



(9.14) 



when Re(l 



Property 9.7. The following list shows how it is possible to express some math- 
ematical functions in terms of the H -function: 



cosz = 



smz 



-'^o.i 



1 rrl,0 
tin 9 



(0,1) 



Vz 



1.0 



(0,l),(i 1) 



(0,1), ("ii) 



2 > 



hi(l + z) 



i^,(z) 



(1,1), (1,1)' 

(1,1), (0,1) 



\z\ < 1 



'0,2 



Vz (9.15) 



^ (^,1),(^,1). 
where the last one is the modified Bessel function of the third kind. 

Property 9.8. The hyper geometric function, when p<qorp = q-\-l with 
< |z| < 1, can always be expressed in terms of the H -function: 



pt^q ■ • ■ , Ol, . . . , Oq, Z) — -^77^ , n„ 



(l-a^,l)? 



(0,1),(1-6„1)? 

(9.16) 

Refer to (0) 

for a more exhaustive hst of relationships between the H- 
function and some other mathematical functions. With ai = Pj = 1, i = 
1, . . . ,p]j = 1, . . . , g all the above properties for i?-functions except Propertv l9.6l 
yield corresponding properties for G-functions due to the relationship (19.51) . 
Before presenting two results from Prudnikov, Brychkov and Marichev (|45l 
P. 346, 368) concerning calculations of integrals involving elementary functions 
and one or two G-functions that are used in this article, we define some more 
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notation as follows. 

m, n,p, g, u, w, w, a; = 0, 1, 2, . . . ; r, s = 1, 2, 3, . . . (r, s are coprime); 

0<m<(7, 0<7i<p, Q < u < X, Q < V < w; 

V* , + ^ , , P + 1 
= u + V , c = m + n : 



(9.17) 



q p 



p-q 



+ 1, p^^d,-^, 



w — a; 



a a + 1 a + £ — I 



For any integer £ > 0, define 
A(£,a): 
and 

A/fi ( \P\ - ( 2± ?2l Hz °1 + -'^ «2 + 1 ftp + 1 

£ ' T 



(9.18) 



, . . . , 



Theorem 9.1 (Prudnikov, Brychkovand Marichev (|45l)). Define notation in I19.17\ I. 
If one of the 35 conditions stated in fjdi , P. 346) holds^ there holds the relation, 



ic^)r 



i /r 



„M „p+a{x-w} — l —a 



'rm+sv,sn-\-su 



dz 



r„r{p~q) 



(9.|20) 



where 



rp-[-sx 



(A(r, ai), . . . , A(r, a„), A(s, 1 - a - di), . . . , 

A(s, I- a - dx), A(r, a„+i), . . . , A(r, Cp)) 



and 



(^^.).,+.» ^ (A(r,6i),...,A(r,6„),A(s,l-a-ci),..., 

A(s, f - a - Ct„), A(r, 6,„+i), . . . , A(r, 6g)). 

Corollary 9.1 (Prudnikov, Brychkov and Marichev (jH)). When a = y^^ , 
w,x,u=l,v = di= 0, ci — (3 in Theorem \9.1[ \9.2U\) reduces to 



f 

Jo 





(«.)?\ 


dz 







G 



A(s,l-«),A(r, (a,)?) 
A(r,(fe,)?),A(s,l~a-/3)^ 
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